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Chapter  I 


INTRODUCTION 

Manufacturing  industries  which  utilize  metal  machining  have  incor¬ 
porated  an  Increasing  level  of  automation  in  recent  years.  Both  direct 
numerical  control  and  computer  numerical  control  systems  can  be  found 
in  many  diverse  areas  of  manufacturing  [2,4].  It  can  be  stated  that 
the  predominant  reason  for  the  use  of  these  automated  systems  has  been 
for  control  of  the  machining  process  and  not  as  a  means  of  integrating 
the  total  system  in  an  optimum  fashion.  A  major  reason  for  this  dis¬ 
crepancy  is  the  lack  of  mathematical  models  which  effectively  optimize 
the  manufacturing  system.  A  consequence  of  this  narrow  use  of  the  com¬ 
puter  in  manufacturing  systems  is  considerable  unused  capacity  which 
could  be  used  for  an  aid  in  increasing  the  productivity  and  cost  effec¬ 
tiveness  of  these  systems. 

An  area  in  which  automation  can  contribute  significantly  to  im¬ 
proving  machining  systems'  effectiveness  is  in  the  area  or  machine  re¬ 
quirements  planning.  In  general,  the  machine  requirements  (MR)  planning 
problem  can  be  defined  as  the  specification  of  the  number  of  each  type 
of  machine  required  in  a  production  process,  or  group  of  processes,  in 
each  period  during  some  planning  horizon  [18].  The  significance  of  this 
problem  Is  highly  dependent  on  the  type  of  production  process  and  its 
complexity.  Previous  research  on  this  problem  has  mostly  been  narrow  in 
scope.  The  production  systems  considered  have  been  relatively  simple 
in  that  little  consideration  is  given  to  the  Interrelationships 
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of  system  functions  and  components.  Attempts  to  realistically  portray 
the  production  system  have  led  to  complex  solution  procedures  such  as 
those  of  references  [6],  [7],  and  [19].  The  nature  of  the  MR  problem 
is  such  that  It  must  be  modelled  as  a  mixed  integer  problem.  Only  in 
this  way  can  the  number  of  machines  decision  be  meaningful.  Integer 
linear  programming  procedures  require  a  generous  amount  of  time  and 
core  for  even  simple  problems.  Therefore,  the  solution  procedures  them¬ 
selves  have  restricted  the  degree  of  complexity  of  the  models. 

The  information  required  to  develop  a  machine  requirements  model 
includes  the  operating  rates  of  the  Individual  machines  within  machine 
centers.  The  cost  associated  with  a  production  rate  and  the  percentage 
of  defects  related  to  production  are  also  data  which  allow  the  formula¬ 
tion  of  a  realistic  model  for  minimizing  cost  in  a  production  system. 
However,  the  operating  rates  are  the  key  element  which  effect  a  viable 
model.  Manufacturing  machinery  is  either  designed  on  an  alternating 
current  (AC)  or  a  direct  current  (DC)  basis.  The  use  of  alternating 
current  usually  results  in  discrete  operating  rates  for  the  machine 
while  direct  current  leads  to  rates  from  a  continuum  within  a  feasible 
range.  The  possibility  of  choosing  from  an  infinite  number  of  rates 
causes  significant  complexity  in  an  optimization  procedure,  particularly 
one  involving  multistage  operation.  It  would,  therefore,  be  a  distinct 
advantage  if  operating  rate  information  for  a  machine  could  be  dis¬ 
cretized  to  further  the  feasibility  of  optimization. 

When  machining  a  particular  product,  the  operating  specifications 
of  a  machine  such  as  the  feed  rate,  depth  of  cut  and  speed  of  the  cutter 
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must  be  optimally  chosen  to  manufacture  a  product  in  a  cost  effective 
manner.  These  operating  parameters  are  what  determine  the  production 
rate  for  a  machine,  the  same  rates  needed  for  the  machine  requirements 
model.  Computer  technology  makes  it  possible  to  more  effectively  deter¬ 
mine  these  specifications,  particularly  where  DC  machines  are  concerned. 
This  class  of  problems  can  be  referred  to  as  machine  parameter  (MP) 
optimization  problems.  Considerable  work  is  documented  in  the  litera¬ 
ture  on  developing  MP  models  but  the  bulk  of  this  research  has  led  to 
models  which  are  nonlinear.  Solution  procedures  such  as  geometric  pro- 
granming  [25]  or  penalty  function  techniques  [30]  have  proven  to  be 
either  extremely  difficult  to  apply  or  are  slow  to  converge  to  an  optimum. 

While  AC  machinery  has  discrete  rates  and  associated  discrete 
machine  parameter  specifications,  an  MR  model  of  broad  applicability  must 
be  able  to  incorporate  DC  machines  as  well.  A  critical  aspect  of  the 
MR  model  is  that  optimizing  each  stage  of  the  process  independently 
does  not  necessarily  optimize  the  entire  system.  The  MR  model,  there¬ 
fore,  needs  complete  operating  rate  information  for  each  machining 
center  (or  stage)  if  the  entire  system  is  to  be  optimized.  Consequently, 
an  associated  area  of  research  for  this  thesis  is  to  develop  a  mathe¬ 
matical  model  for  optimizing  machining  parameters  which  provides  com¬ 
plete  discretized  data.  A  basic  assumption  in  this  respect  of  the 
research  is  that  a  tool  will  only  be  changed  between  passes  of  a  tool 
over  a  workpiece.  This  assumption  reduces  much  of  the  nonlinearity  in 
the  problem  and  more  importantly,  leads  to  a  discretized  solution  pro¬ 
cedure  using  an  efficient  search  algorithm.  With  discrete  operating 
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rates  avail  able,  the  machine  requirements  model  can  be  applied  to  any 
type  of  machining  system,  whether  powered  by  a  discrete  or  continuous 
source. 

The  major  purpose  of  this  research  is  to  develop  a  mathematical 
model  of  the  machine  requirements  problem  which  achieves  the  following 
objectives: 

1.  Describes  an  integrated,  serial  manufacturing  system. 

2.  Can  be  utilized  to  optimize  any  type  of  machining  system, 
whether  of  a  discrete  or  continuous  operating  mode. 

In  addition,  there  is  the  associated  objective  of  developing  a  solution 
algorithm  which  is  efficient  in  both  solution  time  and  core  require¬ 
ments  as  well  as  providing  sensitivity  information  with  respect  to 
final  demand  requirements. 

The  research  will  be  presented  in  a  sequence  that  begins  with  a 
survey  of  the  literature  in  the  areas  of  machine  parameter  optimiza¬ 
tion  and  machine  requirements  planning.  Following  this  is  a  descrip¬ 
tion  of  the  machine  parameters  model  and  solution  procedures.  The 
machine  requirements  model  and  its  solution  procedure  are  then  presented. 
This  formulation  is  tested  and  the  results  are  analyzed  in  relation  to 
a  more  complex  integer  programming  procedure.  Algorithm  descriptions 
and  information  flow  charts  are  provided  to  assist  in  the  understanding 
of  the  models  and  associated  solution  procedures. 


Chapter  II 


LITERATURE  SURVEY 


2.1  Introduction 

This  literature  survey  will  be  presented  in  two  distinct  parts. 

The  development  of  machine  parameter  optimization  techniques  will  be 
discussed  initially.  Then  the  status  of  machine  requirements  models 
will  be  reviewed. 

2.2  Machine  Parameter  Optimization 

There  have  been  many  contributions  to  the  machine  requirements 
optimization  literature  throughout  the  world.  The  most  notable  research 
has  been  carried  out  in  Russia,  Japan  and  the  United  States.  An  obser¬ 
vation  which  should  be  made  at  this  time  is  that  machine  parameter 
models  are  generally  designed  using  one  of  three  criteria:  minimum 
production  cost  per  part,  minimum  production  time  per  part  or  maximum 
profit  rate.  This  review  will  not  concentrate  on  any  particular  cri¬ 
terion  since  any  of  the  criteria  can  be  converted  to  another  with  minor 
modification  of  the  model. 

All  machine  parameter  optimization  modelling  requires  some  simpli¬ 
fying  assumptions.  One  of  the  most  basic  models  was  developed  by  Wu 
and  Ermer  [31]  who  found  the  optimum  machining  parameters  for  a  turning 
operation  without  considering  any  constraints  and  by  varying  one  vari¬ 
able  at  a  time.  Ermer  [9]  solved  the  constrained  machining  problem  in 
an  economic  context  using  geometric  progranml ng .  Two  years  later, 
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Petropoulos  [25]  developed  his  model  for  the  optimal  selection  of 
machining  parameters  using  geometric  programming.  Hati  and  Rao  [12] 
developed  mathematical  programming  models  using  SUMT  techniques  to 
determine  the  optimum  machining  parameters  using  all  three  of  the  cri¬ 
teria  mentioned  above.  Their  research  concerns  deterministic  and 
probabilistic  approaches  and  involves  sensitivity  analysis  in  relation 
to  the  cost  coefficients.  Hati  and  Rao  document  significant  conclusions 
but  their  mathematical  programing  technique,  although  general,  can  be 
complex  to  apply  and  inefficient  in  operation. 

Ermer  [10]  developed  a  Bayesian  approach  to  machining  economics 
by  adaptive  control.  However,  this  approach  requires  precise  monitoring 
of  tool  wear,  a  technique  which  has  had  only  limited  development.  A 
significant  model  on  adaptive  searching  systems  for  machining  was  de¬ 
veloped  by  Lishchinski  and  Moshkov  [17].  Their  model  explores  the 
relationship  between  tool  life,  operating  parameters  and  the  power 
equation  to  develop  a  computerized  search  algorithm  used  to  vary  the 
feed  and  speed  parameters.  Although  appearing  to  be  a  valid  approach, 
this  model  is  highly  complex  and  involves  searching  over  a  wide  range 
of  values  for  the  parameters.  In  addition,  the  algorithm  requires  the 
precise  monitoring  of  tool  wear.  Davis,  Wysk  and  Agee  [8]  developed  a 
simplified  approach  to  determining  optimum  machining  parameters  in  an 
adaptive  control  context.  Their  procedure  used  an  efficient  iterative 
search  algorithm  to  reach  an  optimum,  but  discrete  operating  data  can¬ 
not  be  derived  from  their  model. 

An  automated  planning  process  for  the  optimization  of  machining 
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processes  was  developed  by  Challa  and  Berra  [4].  Their  model  used 
transformed  functions  to  achieve  greater  linearity  and  is  solved  using 
a  gradient  projection  method.  The  well  behaved  nature  of  the  objec¬ 
tive  function  leads  to  an  efficient  solution  procedure  but  the  proce¬ 
dure  cannot  be  extended  to  other  machine  operations  in  such  an  effici¬ 
ent  fashion.  The  model  is  developed  only  for  a  milling  operation. 

Hitomi  [12]  discussed  the  optimal  machining  conditions  for  a  multistage 
system  using  all  three  optimization  criteria.  However,  the  assumptions 
in  his  analysis  were  so  restrictive  as  to  render  the  staging  inconse¬ 
quential.  A  purely  probabilistic  approach  was  taken  by  Iwata,  Murotsu, 
Iwatsubo  and  Fujii  [13]  in  which  chance-constrained  programming  was  used. 
Although  it  was  shown  that  the  optimum  cutting  conditions  are  signifi¬ 
cantly  affected  by  probabilistic  considerations,  the  nature  of  parameter 
estimates  is  inherently  imprecise  and  the  assumption  of  normality  is 
extremely  important  but  substantiated  by  little  empirical  data. 

Several  other  models  have  been  developed  which  have  very  restric¬ 
tive  applications  or  extremely  complex  solution  procedures.  Ravignani 
[26]  attempts  to  solve  optimum  machining  conditions  minimizing  cost 
while  producing  to  a  fixed  demand.  This  is  actually  a  segmented  ap¬ 
proach  where  subsystems  are  isolated  for  individual  optimization  and 
then  the  effects  of  Idle  time  and  other  disturbing  events  are  dis¬ 
cussed.  There  is  no  integrated  systems  approach  for  optimization. 
Boothroyd  and  Rusek  [3]  discussed  the  maximum  rate  of  profit  criteria 
for  determining  cutting  speed.  Their  treatment  results  in  the  conclu¬ 
sion  that  the  maximum  rate  of  profit  criteria  is  a  compromise  between 
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the  conditions  of  minimum  cost  and  minimum  production  time.  The  conclu¬ 
sion  is  also  drawn  that  the  usefulness  of  their  procedure  is  limited  by 
inaccurate  tool  life  data  and  little  improvement  can  be  expected  in 
machining  economics  until  greater  sophistication  is  achieved  in  tool  life 
models.  Claycombe  and  Sullivan  [5]  contribute  the  use  of  response  sur¬ 
face  methodology  to  selecting  a  cutting  tool  in  order  to  maximize  profit. 
This  method  requires  extensive  testing,  data  collection  and  regression 
analysis  and  would  only  prove  feasible  to  companies  with  extensive  use 
of  the  same  cutting  tools  and  the  capability  for  analysis. 

A  comparison  of  nonlinear  programming  techniques  as  applied  to  the 
machine  parameter  optimization  problem  is  provided  by  Kimbler,  Wysk  and 
Davis  [16].  From  this  analysis,  SUMT  and  a  general  penalty  function 
approach  posed  by  the  authors  appear  to  be  the  most  useful;  but  both  are 
limited  by  complexities  in  the  objective  and  constraint  functions  and 
the  choice  of  penalty  term  values.  If  npnl inear  programing  techniques 
are  to  be  more  generally  applicable  to  the  machine  parameter  problem, 
much  additional  research  is  necessary  to  improve  efficiency  and  broaden 
their  feasibility.  Table  2.1  is  a  compilation  of  the  most  significant 
contributions  to  the  MP  problem. 

2.3  The  Machine  Requirements  Problem 

Only  a  limited  number  of  procedures  for  analyzing  the  machine  re¬ 
quirements  problem  have  appeared  in  the  literature.  Most  of  the  studies 
have  been  deterministic,  static  analyses  which  use  a  descriptive  rather 
than  a  normative  approach.  Shubin  and  Madeheim  [29]  developed  a  model 
for  a  single  product,  single  work  center  with  one  operation.  Such  a 
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single  period  model  can  meet  a  prescribed  production  requirement  but 
has  a  very  limited  systems  scope.  Similar  descriptive,  deterministic 
approaches  have  been  reported  by  Apple  [1],  Muther  [23],  Johnson  [15], 
Reed  [27],  and  Moore  [20].  Although  most  of  the  early  work  was  con¬ 
cerned  with  the  single  product,  single  work  center,  single  operation 
situation,  Apple  [1]  and  Francis  and  White  [11]  did  extend  the  analysis 
to  multiple  work  centers.  A  serial,  single  product,  fixed  routing 
operation  is  theorized  but  no  production  system  component  interaction 
is  considered. 

Morris  [21]  developed  a  probabilistic  approach  to  the  MR  problem 
of  a  single  product.  Production  requirements,  operation  performance 
time  and  machine  effectiveness  were  modelled  as  random  variables  but  no 
explicit  form  of  the  distributions  was  developed.  Reed  [28]  also 
treated  the  static,  descriptive  problem  from  a  probabilistic  viewpoint. 
His  multiple  work  center,  serial  flow  system  has  the  objective  of  find¬ 
ing  the  overall  number  of  machines  to  meet  a  production  requirement. 
Specific  distribution  functions  for  the  basic  parameters  were  assumed 
and  not  developed. 

Morris  [22]  contributed  the  earliest  normative  approach  reported 
in  the  literature.  His  decision  model  was  based  on  a  linear  cost  cri¬ 
terion  and  no  constraints.  The  machine  requirements  were  modelled  as 
a  probability  distribution  which  was  unspecified.  All  of  the  models 
to  this  point  have  resulted  in  continuous  solutions  where  the  actual 
number  of  machines  required  results  from  a  round  up  or  round  down  pro¬ 
cedure.  Davis  and  Miller  [6]  approached  the  MR  problem  using 
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deterministic,  steady-state  analysis  in  a  normative  model.  The  formu¬ 
lation  is  based  on  meeting  production  requirements  in  a  single  planning 
period  at  minimum  cost  with  the  number  of  machines  constrained  to  be  an 
integer.  An  Integer  linear  programming  model  was  developed  but  a  deci¬ 
sion  criteria  was  formulated  which  precluded  the  use  of  the  more  complex 
and  costly  ILP  model.  However,  the  decision  tree  type  of  analysis  can 
be  very  tedious  for  multistage  problems.  The  model,  however,  does  in¬ 
clude  tradeoffs  for  investment  cost  and  straight  time  and  overtime 
operating  and  maintenance  costs.  Thus,  MR  models  now  involved  more 
realistic  considerations  and  became  more  feasible  with  the  stagewise, 
Integer  constrained  format. 

Miller  and  Davis  [19]  also  formulated  a  dynamic  resource  alloca¬ 
tion  model  for  the  MR  problem.  Limitations  on  such  resources  as  capital 
budget,  floor  space  and  overtime  hours  were  explicitly  considered  and 
the  interaction  between  machine  centers  was  represented  by  consumption 
of  these  resources.  The  problem  then  becomes  a  deterministic  resource 
allocation  problem  constrained  to  attain  integer  machine  requirements. 

A  drawback  in  the  formulation  is  that  the  results  represent  changes  in 
machine  levels  and  the  initial  number  required  must  be  calculated  to 
satisfy  time  phased  production  requirements.  The  model  is  therefore 
restricted  to  consider  only  a  constant  machine  rate  at  each  center. 

A  significant  step  forward  in  MR  analysis  was  made  by  Davis  and 
Miller  [7].  This  model  formulation  resolves  the  integer  number  of 
machines  required  and  their  discrete  rates  of  operation  for  a  discretely 
distributed  demand  in  a  minimum  cost  framework.  The  model  developed  is 
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for  a  serial,  multistage  manufacturing  system  which  has  both  straight 
time  and  overtime  operating  periods.  Also  provided  for  is  the  deter¬ 
mination  of  in-process  inventories  between  stages.  No  previous  models 
were  designed  to  consider  the  cost  of  variable  operating  rates  and 
integration  of  the  production  and  inventory  systems.  The  model  can  be 
reoptimized  if  the  demand  distribution  changes  in  order  to  update  a 
decision  policy  for  addition  or  disposal  of  machines.  Table  2.2  pre¬ 
sents  a  consolidation  of  major  contributions  to  the  solution  of  machine 
requirements  problems. 

Although  considerable  research  has  been  documented  on  both  the 
machine  parameter  and  machine  requirements  problems,  no  definitive 
analysis  appears  in  the  literature  for  integration  of  the  two  problems 
into  a  single  optimization  model.  The  optimum  cost  of  operating  a 
machining  system  is  directly  related  to  the  minimum  cost  framework  of 
the  machine  parameters  problem.  However,  to  incorporate  the  two  pro¬ 
blems  into  a  single  model  in  the  most  widely  applicable  manner,  requires 
that  both  discrete  and  continuous  machining  systems  be  modelled.  As 
noted  in  the  survey  of  the  MP  literature,  the  nonlinearity  of  the  models 
greatly  complicates  the  solution  procedures  and  a  method  of  discretizing 
all  MP  problems  is  required  if  a  valid  machine  requirements  model  is  to 
be  formulated.  To  achieve  such  an  integration  of  the  two  optimization 
models  would  be  a  major  step  forward  in  optimization  of  the  total  manu¬ 
facturing  system. 
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Chapter  III 


MACHINE  PARAMETERS  OPTIMIZATION  MODEL 


3.1  Introduction 

For  a  machine  requirements  model  to  be  generally  applicable  to  any 
type  of  machining  system,  whether  AC  or  DC,  it  is  necessary  to  determine 
a  set  of  operating  rates,  associated  production  costs  and  the  percent 
of  defects  for  a  particular  operating  rate.  A  machine  parameters  model 
is  to  be  developed  which  achieves  this  information  in  a  form  which  can 
be  readily  applied  to  the  machine  requirements  model. 

The  percent  of  defects  is  analogous  to  the  probability  of  a  defec¬ 
tive  part  and  little  definitive  work  has  been  done  in  this  area.  The 
data  necessary  for  inclusion  of  defectives  in  the  model  must  be  deter¬ 
mined  from:  operational  data  from  a  similar  production  system,  experi¬ 
mental  data  from  sample  production  runs  using  statistical  methods,  or 
incorporated  into  the  model  after  the  actual  production  system  is 
operating.  In  any  event,  the  probability  of  defects  occurring  would  be 
a  highly  dynamic  term  and  would  need  to  be  updated  at  frequent  inter¬ 
vals  to  maintain  viability  of  the  model.  The  cost  of  defectives  will 
be  Initially  defined  in  the  MP  model  formulation  but  will  not  be  pur¬ 
sued  further  In  the  research  except  as  it  relates  to  the  machine  re¬ 
quirements  model.  Table  3.1  presents  a  list  of  notations  used  in 
the  formulation  of  the  machining  parameters  optimization  model. 
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Table  3.1  Notation  for  Machine  Parameters  Optimization  Model 

dc  *  maximum  depth  of  cut  to  be  removed  (mm) 

d  *  depth  of  cut  per  tool  pass  (nro)  where  de{d^ ,  dg#  .  .  .  .  d^} 

f  =  feed  rate  (ran/ rev) 

V  *  cutting  speed  (nm/min) 

L  a  length  of  cut  (mm) 

D  *  diameter  of  cutting  tool  (mm) 

M  =  torque  produced  on  the  tool  (kgM) 

Th  *  tool  thrust  (kg) 

Cj  =  idle  time  cost  ($/min) 

n  =  number  of  passes  made  with  the  tool 

C  *  Taylor  tool  life  constant 

Tr  *  time  required  to  replace  a  tool  (min) 

Tl  *  tool  life  (min) 

P  =  power  at  the  tool  (kw) 

Cp  *  cost  per  part  ($) 

Cq  *  cost  of  cutting  ($/min) 

Cj  *  cost  of  changing  a  tool  ($/min) 
np  *  number  of  tool  passes  per  part 
nt  *  number  of  tool  changes  per  pass 
Tp  *  time  per  pass  (min) 
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3.2  Problem  Specification 

Specifically,  the  machine  parameter  optimization  problem  is  one 

<• 

of  determining  the  operating  parameters  of  feed,  speed,  and.dep'th  of  cut 
in  order  to  optimize  a  particular  manufacturing  objective.  The  objec¬ 
tive  considered  in  this  research  is  one  of  minimum  cost  of  production 
per  part.  The  total  cost  to  machine  one  part  (Cp)  can  be  represented 

*  v 

by  the  following  expression, 

CP  =  CITI  +  CcTc(V*f’d)  +  cTTT(V»f*d)  +  c/d(V’f’d)  (3J) 

where , 

Tj  =  Machine  idle  time  due  to  setup,. load  and  unload  operations  and 
idle  tool  motions,  a  constant. 

Tc  =  Time  In  cut  which  is  the  expected  time  per  pass  times  the  re¬ 
quired  number  of  passes. 

Ty  a  Tool  change  time  which  is  the  expected  number  of  tool  changes 
per  part  times  the  replacement  time. 

Pd  a  Probability  of  a  defective  part  as  a  function  of  V,  f  and  d. 
Cj,  Cc,  CT,  Cd  =  Cost  per  minute  for  idle  time,  cutting  time,  tool 
replacement  time  and  defectives,  respectively. 

The  number  of  passes  per  part  is  restricted  to  being  an  integer 
and  is  defined  by  the  expression, 
d 

np  a  -g-,  an  integer  where  de{d^ ,  d^,  .  .  .  ,  dj^ > .  (3.2) 

Based  on  the  assumption  that  it  is  Infeasible  to  change  a  tool  during  a 
pass,  the  number  of  tool  changes  per  pass  can  be  defined  as, 


„  _  time  per  pass  _  P 
nt  tool  life  Tj^ 


(3.3) 


where  both  numerator  and  denominator  are  bounded  and  the  feasible  set 
of  n^  values  is  lower  bounded  by  an  integer  fractional  value  and  upper 
bounded  by  one.  Thus, 


nt  =1,  (3.5) 

max 

since  tool  life  is  required  to  be  at  least  equal  to  the  time  required 
for  one  pass.  The  number  of  tool  changes  per  part  is  defined  by, 

nT  *  np  •  nt  (3.6) 

and  Ty  =  nT  •  Tr,  where  TR  is  the  tool  replacement  time,  or 

tt  =  V  nt  *  tr  (3-7) 

Utilizing  these  relationships,  the  objective  function  (3.1)  can  now  be 
written  as, 

Cp  *  ^I7I  +  *  7P  *  np  +  ^T  ’^R  *  nt  ’  np  (3.8) 

subject  to  a  set  of  constraints  which  will  be  illustrated  here  for  turn 
ing  and  boring  of  mild  steel  alloys.  The  origin  of  these  constraints 
will  not  be  discussed  here  but  adequate  development  appears  in  other 
sources,  such  as  references  [1]  and  [10].  The  constraining  relation¬ 
ships  are: 
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n  e  {Np},  where  Np  is  the  feasible  set  of  n  values, 

H  r 

e  {Nt>  •  where  Nt  is  the  feasible  set  of  nt  values. 


^min  —  ^  —  *max 


Vm1n  i  V  i  ^max 


T.  <  TL  <,  T. 
Lmin  —  Lmax 


M  =  0.084  f*75  •  d  *  D  <  M 

—  max 


Th  =  2.02  f JZ  •  d  £  Th 


max 


P  =  0.746  M  •  V/4500  <  P 


max 


SF  =  g(V,f.d)  <  SF, 


max 


limit  on  spindle  speed 
limit  on  feasible 

limit  on  spindle  torque 
limit  on  tool  thrust 

limit  on  available  power 

specific  functional  limit 
on  achieving  a  desired 
surface  finish 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


The  optimization  model  will  consider  only  initial  cuts  up  to  the 
final  pass;  and,  therefore,  surface  finish  considerations  will  not  be 
included.  With  these  considerations  and  the  realization  that  idle  time 
is  a  constant,  the  objective  function  can  now  be  condensed  to, 


Cp  =  (np  •  o^)(Cq  ■  +  Cj  •  T^) 

Utilizing  the  familiar  expanded  Taylor  tool  life  equation, 

T  *  - £ - 

L  Va  •  f0  •  dY 


(3.16) 


(3.17) 


and  the  time  required  for  a  single  pass  expression, 
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T  _  IT  •  D  •  L 

TP  “V"."7  * 

the  following  relationships  can  be  developed: 


(3.18) 


V  = 


C1/a 

T^l/a  .  f&/a  .  dy/a  ’ 


(3.19) 


and,  TL  3  ’p  a  I  y  ,  or  in  expanded  form  as 

L  t 


tl  = 


tt  •  D  •  L  •  nt_1  •  TL1/a  •  fB/a'}  •  dY/a 
- ? 


,1/a 


(3.20) 


Solving  for  feed  (f),  we  obtain, 


f  * 


Hr  (1-1/a)  . 

J. _ 

T  •  U  • 


r/alcT57^TTJ 


(3.21) 


Consequently,  for  any  given  value  of  np,  n^.,  and  T^,  the  machining  feed 
and  speed  can  be  determined  from  an  optimization  procedure  developed 
using  these  relationships.  The  objective  function  (3.16)  is  to  be  mini¬ 
mized  subject  to  a  set  of  constraints  given  by  equations  (3.9)  through 
(3.14)  and  the  integer  restrictions  of  (3.2)  and  (3.3).  The  model  is  now 
In  a  form  whereby  an  iterative  search  procedure  can  be  used  to  take  ad¬ 
vantage  of  the  discrete  nature  of  the  variables  nt  and  np  and  the  sec¬ 
tional  convexity  of  the  objective  function. 
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3.3  Illustrative  Example 

The  machine  parameter  optimization  problem,  as  defined,  has  been 
solved  using  an  example  problem  with  the  following  data: 

D  =  100  mm,  L  =  1000  mm,  dc  =  5  mm,  d^  =  1.2  mm,  dmax  *  2.75  mm, 

fm1n  "  *3  mm/rev»  fmax  *  *75  mn/rev’  Vmin  “  5  m/min» 

V  *  4  x  105  mm/min,  Pmav  =  2.25  kW,  C  =  6  x  1026,  T.  *25  min, 

max  max  Lmin 

T,  *  45  min,  a  =  5.0,  6  =  1.75,  y  =  0.75,  TD  =  1.5  min, 
max  K 


T.  *  60  kg,  Mmav  =  47.5  kgM,  Cr  *  .3  $/min,  CT  *  .5  $/min. 
nmax  max  u  1 


The  basis  for  this  example  data  is  a  carbide  tool  and  mild  steel  work- 
piece  upon  which  a  turning  operation  is  performed.  In  observing  the 
problem  constraints,  it  can  be  seen  that  substitutions  can  be  readily 
made  using  equations  (3.2),  (3.17),  and  (3.21).  These  substitutions  re¬ 
sult  in  constraint  equations  written  in  terms  of  nt,  np,  and  T^.  Since 
nt  and  np  are  integer  terms,  the  most  simplifying  step  is  to  rewrite 
all  of  the  constraints  in  terms  of  TL  as  the  bounded  variable.  With 
constraints  (3.9),  (3.10),  (3.12),  (3.13),  and  (3.14)  altered  in  this 
way  and  (3.2),  (3.3)  and  (3.11)  left  unaltered,  the  problem  can  now  be 
written  in  a  multistage  decision  context  with  an  explicit  optimization 
rationale  at  only  one  stage.  This  staged  representation  is  illustrated 
in  Figure  3.1. 

In  this  staged  context,  S^<*>  represents  the  vector  of  bounding 
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limits  on  tool  life  given  by  the  altered  constraints.  Sg<np>  represents 
the  same  vector  altered  by  the  assignment  of  np  from  its  feasible  set 
and  likewise  for  S.i<nt>  and  n^.  Consequently,  S^<nt>  is  a  vector  of 
upper  and  lower  bound  values  for  T^,  of  which  max  {S^}  represents  the 
maximum  lower  limit  attainable  by  and  min  {S^}  represents  the  cor¬ 
responding  upper  bound  limiting  value. 

A  further  limiting  of  the  problem  is  attained  by  establishing  a 
more  restrictive  lower  and  upper  bound  on  nt  through  analysis  of  the  tool 
life  equation  and  the  defining  relationship  for  nt>  Since  the  tool 
life  ratio  and,  therefore  the  nt  ratio.  Is  most  sensitive  to  V,  the 
minimum  attainable  V  and  nt  can  be  derived  in  the  following  manner: 


min 


f8 

l  wax 


dl 


'max 


max 


1/a 


(3.22) 


nt  =  TP/Tl  «  ,  P-*-L 

Srin  Lmax  T 


L  ’  vm1n'  *  fmax  [T,  •  V  .  •  f 

max  I  LMax  mm  i 


(3.23) 


max 


it  •  D  •  L 


in  a  similar  manner  n^x  can  be  defined, 

1 


n. 


max 


r  •  v  .  f  3 

Lmin  max  min 


•  6  •  L 


(3.24) 


Figure  3.2  is  a  recapitulation  of  the  mathematical  model 
written  In  terms  of  np,  nt  and  T^.  Figure  3.3  shows  this  model  with 
the  above  data  inserted. 


Hlnlatie  Co  *  (fln  '  nt^c 


Figure  3.2:  Resultent  model  structure 
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Minimize  Cp  *  (np  •  nt)(0.3  Tj_  +  0.75) 
Subject  to: 

np  e  {2,3,4} 

n  ,111  It 

nt  E  T5’’  S’ . T} 


3.8574 

1.25 

•»  T  > 

2.1264 

n  n  *15 

nt  ’  nP 

-  L  - 

l  rnr 

"t  •  "p 

1.25 


0.4454  •  n, 


.15 


-1-0.97 


(0.567  •  nt  •  n  *15) 


TO 


35.6323  *  n. 


.15 


(0.567  •  nt  •  n  *15) 


-T75T 


25  <  \  <  45 


TL^ 


(5.9406  •  n  )”'90  (-1  v7^3.7— ) 

p  "t'V  . 


1.25 


TLi 


/•l  i _  \-.867  /  1  .7637  i 

d.131  •  n  j  ( - 


nt  ’  nP 


1.25 


TLi 


(551.29  •  np",85)(0.567  •  nt  •  np-'3) 


.15 


.4 


-8.33 


Figure  3.3:  Mathematical  model  for  example  problem. 
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An  Iterative  search  algorithm  was  developed  to  take  advantage  of 
these  discrete  aspects  of  the  problem  formulation.  Figure  3.4  is  an 
Information  flow  chart  of  the  algorithm  and  Table  3.2  contains  solution 
information  for  the  example  problem.  The  optimum  solution  is  at 
iteration  4  and  is  as  follows: 

d*  *  2.5,  f*  «  0.75,  V*  =  114.232,  T,  =  25.0000. 

LP 

The  execution  time  required  to  produce  all  feasible  solutions  was  0.27 
sec  using  a  FORTRAN  algorithm  Implemented  on  an  IBM  370/Mod  165  com¬ 
puter.  However,  taking  advantage  of  the  sectional  convexity  of  the 
objective  function,  termination  of  the  solution  procedure  upon  reach¬ 
ing  a  minimum  reduces  the  execution  time  to  0.11  seconds. 

Table  3.3  presents  a  relative  comparison  of  the  solution  results 
obtained  for  the  above  example  problem  when  solved  using  several  exist¬ 
ing  algorithms  and  the  staged  approach  given  here.  Not  only  does  the 
approach  given  here  produce  an  optimum  result  much  faster;  but,  more 
importantly,  it  is  the  only  algorithm  which  generates  a  solution  that 
reflects  tool  changes  only  between  passes  (an  important  practical  con¬ 
sideration)  . 

3.4  Conclusions 

The  realistic  assumptions  utilized  in  developing  the  machine  para' 
meter  optimization  model  presented  here  offer  the  basis  for  a  much  more 
efficient  algorithm  than  previously  available.  The  method  is  very 
amenable  to  planning  optimum  machine  settings  and  can  be  used  to 
initialize  an  adaptive  control  system  which  will  then  reoptimize  feed 


Table  3.2  Feasible  Results  for  Example  Problem 
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Table  3.3  Comparison  of  Solution  Times  for  Example  Problem 


A1 qorl thm/Ref erence 


Execution  Time  (sec.) 


SUWT  [16]  4.88 

Kimbler,  Wysk,  Davis  [16]  2.85 

Davis,  Wysk,  Ages  [8]  0.46 

Hayes,  Davis  (Thesis  Algorithm)  0.11 
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and  speed  during  cut,  as  described  in  [8].  However,  the  principal  ad¬ 
vantages  of  the  model  are:  the  discrete  nature  of  the  resultant  solu¬ 
tion  process,  consequent  reduction  in  solution  time,  and  improved  solu¬ 
tion  accuracy.  The  discrete  parameter  information  of  Table  3.2  can  now 
be  incorporated  into  a  generally  applicable  machine  requirements  model. 
The  production  rate  data  needed  for  the  MR  model  can  be  derived  from 
the  machine  parameters  model  by  calculating  the  production  time  per 
unit  which  is  imbedded  in  the  formulation.  The  operating  rate  for  the 
machine  is  then  the  reciprocal  of  this  production  time. 


Chapter  IV 


DEFINITION  OF  THE  MACHINE  REQUIREMENTS  MODEL 
4  .1  Introduction 

The  problem  addressed  in  this  section  is  one  of  defining  the 
optimum  number  of  machines,  their  operating  rates  and  the  time  of 
operation  at  each  rate  to  accomplish  the  production  of  a  specific 
product.  The  manufacturing  system  in  which  production  occurs  conforms 
to  a  serial,  multistage  configuration  with  a  different  processing 
operation  occurring  at  each  machining  center  (or  stage).  It  is  as¬ 
sumed  that  product  demand  and  production  are  given  on  a  per  day  basis. 
Within  this  general  framework,  a  model  is  developed  which  can  be  em¬ 
ployed  to  answer  the  following  fundamental  questions  related  to  the 
design  of  such  a  manufacturing  system,  in  the  context  of  cost  mini¬ 
mization. 

1.  How  many  machines  are  needed  at  each  machining  center  of 
the  system? 

2.  What  are  the  optimum  operating  rates  for  the  machines  in 
each  production  stage? 

3.  How  long  should  the  machines  be  run  in  each  stage  for  a 
specific  final  production  level? 

Approaches  to  the  first  of  these  questions  were  addressed  in 
Chapter  II  but  the  fundamental  framework  of  the  model  will  closely  fol¬ 
low  that  of  Davis  and  Miller  [7].  Such  a  methodology  treats  the 
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manufacturing  process  as  an  integrated  system  where  production  costs 
vary  with  processing  rates.  The  rates  at  each  stage  exist  as  either  a 
continuum  over  a  range  (for  DC  type  machinery)  or  as  a  discrete  set  of 
feasible  rates  (for  AC  machines).  However,  the  methodology  discussed 
in  Chapter  III  allows  a  discrete  set  of  operating  rates  for  use  in  the 
machine  requirements  model  whether  AC  or  DC  machinery  is  employed. 
Variation  in  production  cost  as  the  rate  varies  is  due  to  such  factors 
as  tool  wear,  power  consumption,  maintenance  cost  and  direct  labor, 
factors  which  were  integral  elements  of  the  machine  parameters  model. 

It  is  assumed  that  the  percent  defectives  resulting  from  a  manufactur¬ 
ing  process  is  a  function  of  the  processing  parameters  as  discussed 
in  Section  3.2.  Therefore,  a  defectives  rate  can  be  associated  with 
each  discrete  processing  rate. 

4.2  Definition  of  the  Mathematical  Model 

Three  basic  assumptions  underlie  the  machine  requirements  model: 

1.  The  machines  in  each  work  center  are  homogeneous. 

2.  Each  production  day  is  broken  down  into  two  components,  an 
8  hour  first  shift  and  an  8  hour  second  shift  period. 

3.  The  output  from  a  production  stage  dictates  the  amount  to  be 
processed  at  the  next  stage. 

A  list  of  notations  used  in  the  development  of  the  machine  requirements 
model  is  contained  in  Table  4.1. 

For  the  above  system  description,  the  objective  for  a  minimum 
cost  configuration  is  one  which  minimizes  the  following  daily  operat¬ 
ing  cost  function: 
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Table  4.1  Notation  for  the  Machine  Requirements  Model 

N  =  the  number  of  machine  centers  in  the  system 

r.j  =  a  feasible  operating  rate  from  a  set  R  at  stage  i 

C*  =  manufacturing  cost  for  one  hour  of  production  at  rate  ri 
during  the  first  shift  ($/hour) 

C°  *  manufacturing  cost  for  one  hour  production  at  rate  r-  during 
i  1 

the  second  shift  ($/hour) 

br  *  percent  of  defective  units  incurred  at  stage  i  by  processing 
at  rate  r^ 

CF  =  fixed  cost  of  a  machine  per  day  in  stage  i  ($/day) 
i 

n^.  =  number  of  machines  utilized  in  stage  i 

t^  *  number  of  hours  of  operation  at  rate  r.  during  the  first  shift 

t°  =  number  of  hours  of  operation  at  rate  r-  during  the  second  shift 

ri  1 


time  period. 
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N  e  e  A  n 

Total  Cost  *  i  l  (Cj  t*  +  C°  t°  +  CF  n  )  (4.1) 

1=1  reRi  r1  ri  ri  ri  hi  1 

processing  cost  equipment  cost 

subject  to  the  following  restrictions,  for  i  =  1,  .  .  .  ,  N. 

First,  the  total  number  of  units  processed  at  stage  i  must  equal 
the  quantity  available  for  processing  at  that  stage, 

I  r.(t*  +  t°  )  *  S,  (4.2) 

reR,  ’  r1  r1  1 

Second,  the  quantity  of  output  product  at  a  stage,  quantity  pro¬ 
cessed  less  the  fraction  defective,  must  equal  the  quantity  processed 
by  the  next  stage, 


4  r'(1  •  W  *?,’  * 5'  ■ s" 


(4.3) 


For  the  final  stage  this  restriction  takes  the  following  form, 


l  rn(l  -  bf  )(tj  +  t°  )  -  Sn  =  final  demand 


reR 


(4.4) 


i 


Finally,  the  units  being  processed  at  a  stage  cannot  employ  more 
processing  time  than  is  available  on  the  number  of  machines  allocated 
to  that  stage.  For  the  first  shift  operation  this  is: 


t 

j 


I 

1 


l  tl  l8n1  (4.5) 

reRt  ri 


■A*.  • 


I 


and  for  the  second  shift  operation. 


I  tj  <  8n, ,  and  (4.6) 

reR,j  i  1 

.  t°  ,  n.  >  0  and  n.  an  integer  (4.7) 

The  above  model  is  a  fundamental  linear  program  requiring  an 
integer  domain  for  the  subset  of  variables  associated  with  the  number 
of  machines  and  allowing  the  other  decision  variables  to  have  a  con¬ 
tinuous  solution.  A  multistage  system  employing  this  model  can  be 
solved  using  a  mixed  integer  linear  programming  algorithm  with  a  branch 
and  bound  solution  technique.  However,  the  solution  time  and  core  re¬ 
quirements  for  such  a  procedure  can  become  prohibitive  for  a  problem 
involving  a  significant  number  of  stages  and  operating  rates.  An 
alternative  approach,  using  dynamic  programs ng ,  may  offer  distinct 
computational  advantages  and  result  in  a  greater  variety  of  decision 
information  for  the  analyst.  Such  an  approach  will  be  formulated  in 
the  next  chapter. 


Chapter  V 


DYNAMIC  PROGRAMMING  SOLUTION  PROCEDURE 

5.1  Rationale  of  a  Dynamic  Programing  Formulation 

The  staged  character  of  the  manufacturing  system  being  modelled  and 
the  nature  of  the  transition  between  stages  as  described  in  equations 
(4.2)  and  (4.3)  lead  directly  to  investigating  a  solution  procedure  using 
dynamic  programming.  Figure  5.1  depicts  a  multistage  serial  production 
system  with  decisions  at  each  stage  represented  by  the  number  of  machines 
and  operating  hours  decision  variables.  The  return  ,  represents  the 
cost  incurred  at  each  stage  in  producing  to  a  specified  final  demand  SN. 
The  transition  between  stages  is  represented  by  the  state  variable  . 
Solving  for  the  necessary  S  value  at  each  stage  dictates  the  number  re¬ 
quired  as  output  from  the  previous  stage.  Through  working  backwards 
through  the  problem,  a  set  of  decisions  can  be  determined  at  each  stage 
to  attain  the  required  final  production  level  of  S. 

Handling  the  problem  In  the  above  fashion  restricts  the  decisions 
to  those  which  are  tied  directly  to  the  final  demand  value.  However,  a 
much  more  broadly  based  solution  procedure  can  be  achieved  by  using 
state  Inversion  techniques  [24]  thereby  reversing  the  state  transition 
arrows  in  the  serial  system  as  shown  in  Figure  5.2.  The  problem  can 
then  be  treated  as  a  fundamental  dynamic  programming  problem  using  a 
recursive  procedure  to  work  against  the  flow  of  the  arrows.  The  deci¬ 
sion  resulting  from  such  an  inversion  procedure  will  reflect  a  broad 
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range  of  system  states  in  relation  to  both  initial  and  final  values  in 
the  production  process. 

As  the  model  is  basically  a  mixed  integer  linear  program  which  is 
decomposable  into  stages,  it  is  necessary  to  use  linear  programming 
solution  procedures  at  each  stage  of  the  dynamic  program  to  determine 
the  optimum  decisions  which  hold  for  that  stage  and,  therefore,  all  prior 
stages.  Insuring  that  both  optimality  and  feasibility  are  attained  in 
relation  to  the  decisions  at  each  stage  will  assure  that  decisions  for 
the  entire  system  are  optimized  in  accordance  with  the  principle  of 
dynamic  programming. 

5.2  Maintaining  Optimality  and  Feasibility  Through 

Parametric  Analysis 

The  linear  programming  approach  used  in  each  stage  of  the  manufac¬ 
turing  model  Is  actually  one  in  which  the  initial  decision  values  are 
predefined  to  maintain  optimality  and  further  iterations  maintain  opti¬ 
mality  of  the  solution  while  determining  the  range  of  production  values 
for  which  the  current  solution  remains  feasible.  Given  a  set  of  produc¬ 
tion  rates,  unit  cost  information  and  percent  of  defectives  for  the 
Initial  stage  in  the  system,  it  is  possible  to  determine  the  lowest  cost 
production  rate  for  the  first  8  hour  shift.  This  relationship  is  of 
the  form: 


Minimum  cost  rate  =  Min  (C„  /r, (1-b  )) 

r^eRi  rl  1  rl 


(5.1) 


Figure  5.3a  illustrates  the  formation  of  a  linear  programming 


tableau  with  the  three  rows  of  the  constraint  matrix  representing  the 
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restrictions  of  equations  (4.3),  (4.5)  and  (4.6)  of  the  mathematical 
model.  The  right  hand  side  of  the  first  row  of  the  matrix  is  the  pro¬ 
duction  level  for  8  hours  using  one  machine  at  the  lowest  cost  produc¬ 
tion  rate.  The  second  and  third  row  right  hand  side  values  are  the  first 
and  second  shift  slack  time  hours. 

The  tableau  of  Figure  5.3b  illustrates  the  solution  after  an  ini¬ 
tial  pivoting  operation  where  the  variable  associated  with  the  minimum 
cost  rate  enters  the  basis  replacing  the  artificial  variable  A.  The 
production  range  for  which  this  solution  is  feasible  is  zero  up  to  the 
level  determined  by  8  hours  at  r^  .  .  The  key  element  of  linear  pro¬ 
gramming  which  has  repeated  application  with  this  model  is  that  of  post¬ 
optimality  analysis  which  can  be  done  to  the  right  hand  side  in  relation 
to  changes  in  production  level  represented  by  the  first  row.  The  change 
in  production  level  for  which  a  solution  is  feasible  is  determined  by 
the  effect  of  the  first  column  of  the  basis  inverse  on  the  right  hand 
side.  The  feasible  range  of  production  is  then  calculated  by  moving 
this  basis  inverse  column  to  the  right  hand  side  and  solving  for  the 
number  of  units  which  maintain  feasibility.  Each  subsequent  iteration 
through  the  simplex  procedure,  in  which  a  new  variable  enters  the  basis, 
will  maintain  optimality  while  the  artificial  variable  column,  or  first 
column  of  the  basis  inverse,  establishes  the  production  range  for  which 
the  solution  holds.  For  example,  with  the  first  column  of  the  basis 
Inverse  and  the  right  hand  side  vector  having  the  following  values, 
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.0 

8.0 

8.0 

.034 

and 

.0 

leads  to 

.034  b 

-.034 

8.0 

8-. 034  b 

A 

RHS 

RHS  +  Ab 

the  lower  range  for  this  solution  is  the  8  hour  production  rate  deter¬ 
mined  by  the  time  variable  associated  with  the  first  row  and  the  upper 
bound  is  determined  by  the  point  at  which  the  solution  becomes  infeas¬ 
ible.  Solving  the  equation, 

8-0.034  Ab  =  0  (5.2) 

determines  the  Ab  value  or  the  additional  production  above  the  initial 
right  hand  side  first  row  value  for  which  the  current  solution  holds. 
Continuing  to  iterate  while  bringing  in  the  next  most  cost  effective 
production  rate  variables  will  maintain  optimality;  while  determining 
the  effective  range  of  production  through  postoptimality  analysis  of 
the  right  hand  side.  These  iterations  are  accomplished  through  dual 
simplexing  on  the  right  hand  side  value  which  will  first  cause  the  solu¬ 
tion  to  become  infeasible  when  parametric  analysis  Is  used.  Progressing 
In  this  manner  will  Insure  that  there  are  no  discontinuities  In  the 
ranges . 

The  Iterations  at  each  stage  continue  until  either  no  variable 
can  be  brought  into  the  basis  or  the  upper  limit  of  production  is  reached. 
If  no  additional  variable  can  enter  the  basis,  it  is  necessary  to  add  a 
machine  to  accomplish  higher  production  levels.  This  is  done  through 
increasing  the  Initial  hours  available  for  first  and  second  shift 
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production  by  8  for  each  shift  and  continuing  to  iterate  after  insuring 
that  the  solution  is  returned  to  feasibility.  When  the  upper  limit  of 

t 

planned  production  is  attained  at  a  stage,  it  is  necessary  to  continue 
to  the  next  stage. 

5.3  Transition  Between  Stages  of  the  Manufacturing  Process 

Following  the  initial  stage  of  production,  each  stage  of  the  pro¬ 
duction  process,  and  therefore  the  model,  must  absorb  the  returns  (or 
cumulative  costs)  of  stages  previous  to  it.  It  must  also  reflect  the 
state  variable  production  quantity  transition  value  between  stages.  The 
addition  of  a  new  row  to  the  constraint  matrix  representing  equation 
(4.2)  of  the  model  enables  both  of  these  requirements  to  be  met.  Figure 
5.4a  is  the  initial  simplex  tableau  for  all  stages  other  than  the  first. 
S.j  is  the  transition  value  and  represents  the  input  to  stage  i  to  obtain 
the  level  of  production  indicated  for  the  current  solution.  It  is  also 
the  output  needed  from  stage  i-1.  The  cost  of  production  from  previous 
stages  must  be  absorbed  into  the  current  stage.  The  cost  of  production 
at  a  stage  is  represented  by  the  objective  function  value  associated  with 
a  feasible  tableau  and  contains  two  components.  One  is  the  fixed  compo¬ 
nent  consisting  of  the  basic  variable  cost  values  multiplied  by  the  right 
hand  side.  The  fixed  component  also  includes  the  fixed  cost  of  a  machine 
times  the  number  of  machines.  The  variable  component  is  the  basic  vari¬ 
able  cost  vector  times  the  first  column  of  the  basis  Inverse  to  account 
for  the  cost  of  production  above  the  lower  bound  of  the  range  for  which 
the  solution  holds.  In  summary,  this  relationship  can  be  stated  as, 


Figure  5.4b:  Initial  IP  tableau  for  stages  other  than  first  after  first  pivot  operation 
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Total  Cost^  =  Cp  •  n<  +  £Cg][b]  +  [CgltB'1]  (5.3) 

When  a  stage  absorbs  the  cost  from  a  previous  stage,  it  must  absorb 
both  the  fixed  and  variable  components.  The  fixed  portion  is  added 
directly  to  the  objective  function  value  of  the  next  stage  while  the 
variable  component  becomes  the  cost  coefficient  associated  with  the 
transition  function  value  .  Since  Si  is  always  in  the  basis,  the  cost 
coefficient,  C$^,  is  always  included  as  a  component  of  Cg,  the  basis 
cost  vector.  Figure  5.4b  is  the  initial  tableau  of  the  ith  stage  (for 
i  greater  than  1)  following  the  first  pivot  operation,  or  iteration  up¬ 
date.  Thus  the  cost  (or  objective  function  value)  includes  the  cost 
associated  with  the  previous  stage  and  all  stages  prior  to  it.  The  S- 
variable  in  the  basis  gives  the  input  level  needed  at  the  ith  stage  and 
the  necessary  output  from  stage  i-1. 

Since  the  costs  and  transition  values  are  optimized  at  each  stage, 
and  therefore  independent  of  all  previous  stages  after  absorption,  the 
number  of  stages  has  no  effect  on  the  procedure  except  in  total  solution 
time. 

The  most  critical  aspect  of  absorbing  previous  stage  costs  is  that 
of  maintaining  continuity  of  production  ranges.  Cost  functions  which 
hold  for  particular  ranges  in  stage  i  must  be  absorbed  into  stage  i  +1 
with  full  cognizance  of  the  ranges  involved.  A  change  in  cost  function 
in  stage  i  directly  causes  a  change  in  the  fixed  and  variable  portions 
of  the  cost  components  absorbed  into  stage  i  +1.  A  range  split  will 
then  result  in  stage  i  +1.  The  solution  vector  for  stage  i  +1  will  not 
change,  but  the  cost  associated  with  the  solution  changes,  necessitating 
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the  split  in  ranges  to  maintain  correct  accountability  of  absorbed 
costs.  An  element  which  must  be  considered  in  handling  the  absorption 
of  returns  caused  by  a  range  split  is  that  the  range  from  a  previous 
stage  affects  the  current  stage  through  the  transition  function;  and  is 
not  directly  analogous  to  the  current  production  of  good  units  but  to  the 
input  units  to  the  current  stage.  For  instance,  if  a  change  in  returns 
occurred  at  a  certain  production  level  in  the  previous  stage,  that 
change  is  reflected  in  the  current  stage  at  a  production  level  which  is 
reduced  to  account  for  defectives  in  the  production  process. 

5.4  Incorporating  Economic  Tradeoffs  as  Decision  Criteria 

A  key  decision  point  in  the  analysis  of  the  manufacturing  process 
is  whether  it  is  more  economically  advantageous  to  bring  in  higher  cost 
production  rates  or  to  increase  the  number  of  machines  in  order  to  in¬ 
crease  the  level  of  production  in  a  stage  of  the  machining  system.  In¬ 
formation  included  in  the  linear  programming  solution  provides  a  readily 
available  means  of  economic  analysis  to  make  this  decision.  Associated 
with  the  slack  values  in  the  constraint  matrix  of  the  LP  procedure  are 
reduced  cost,  or  "shadow  price,"  values  represented  as  the  Zj-Cj  values 
in  the  cost  row  pertaining  to  the  slack  variables.  These  values  indi¬ 
cate  the  amount  of  savings  which  could  be  incurred  for  each  additional 
hour  of  time  available  in  the  first  or  second  shift  respectively.  It 
must  also  be  considered  that  adding  hours  of  production  can  only  be  done 
in  16  hour  increments  per  machine.  Therefore,  what  appears  as  a  simple 
check  to  determine  If  the  savings  available  exceeds  the  cost  of  adding 
another  machine,  becomes  more  complex  with  the  realization  that  a  full 
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16  hours  of  production  may  not  be  needed  from  another  machine.  Conse¬ 
quently,  the  decision  must  be  based  on  production  levels  as  well  as  costs 
and  savings,  a  task  for  which  the  LP  procedure  is  particularly  well 
suited. 

Through  the  use  of  the  objective  function  cost  values,  it  is  pos¬ 
sible  to  determine  whether  it  is  economically  feasible  to  add  a  machine 
as  well  as  the  precise  production  level  at  which  this  occurs.  Within  a 
particular  stage,  the  reduced  costs  associated  with  the  slack  variables 
are  used  as  an  indicator  of  whether  the  economic  situation  bears  fur¬ 
ther  analysis.  If, 

8  •  (2„  -  CR  )  +  (ZR  -  CR  )  >  C,  ,  (5.4) 

K-j  •'■j  *2  k2  r  i 

it  is  necessary  to  undertake  further  analysis.  The  reduced  costs  are 
multiplied  by  8  hours  since  adding  a  machine  would  provide  two  addi¬ 
tional  shifts  at  8  hours  each  and  the  relationship  of  equation  (5.4) 
is  therefore  balanced. 

Upon  determination  that  it  may  be  advantageous  to  add  another 
machine,  a  formal  analysis  is  conducted  to  determine  if  the  pertinent 
production  level  for  a  breakeven  point  is  within  the  range  associated 
with  the  current  solution  vector.  This  analysis  proceeds  as  follows, 
let  RVCUR  be  the  current  variable  portion  of  the  cost  function  and  RFCUR 
be  the  fixed  portion  as  they  pertain  to  a  basic  LP  solution  with  re¬ 
duced  costs  satisfying  equation  (5.4).  In  order  to  test  whether  adding 
a  machine  is  valid,  it  is  now  necessary  to  hypothetically  do  so  by 
adding  8  hours  to  the  first  and  second  shifts. 
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The  lower  bound  of  the  range  for  which  the  basic  solution  holds  Is 
used  as  the  new  initial  right  hand  side  production  value.  Using  the  dual 
simplex  procedure,  the  feasibility  of  the  solution  can  be  restored  and 
the  lower  bounds  of  the  base  solution  and  the  updated  solution  will  cor¬ 
respond,  although  the  upper  bounds  of  the  ranges  may  differ.  The  objec¬ 
tive  functions  of  the  two  solutions  will  differ  since  they  involve  dif¬ 
ferent  bases  and  this  serves  as  the  means  of  comparison  for  purposes  of 
a  decision. 

Through  setting  the  current  objective  function  of  the  base  solution 
equal  to  the  objective  function  of  the  feasible  solution  with  a  machine 
added,  represented  by  fixed  and  variable  components,  RFP  and  RVP,  as 

RFCUR  +  RVCUR  (Ab)  =  RFP  +  RVP  (Ab)  +  CF  ,  (5.5) 

hi 

a  Ab  value  can  be  determined  for  which  the  cost  functions  are  equal.  This 
indicates  the  production  level  at  which  adding  a  machine  is  economically 
feasible.  If  this  production  level  is  outside  the  range  of  production 
for  which  the  base  and  revised  solution  hold,  then  no  machine  is  added 
and  further  iterations  proceed  from  the  base  solution.  If  adding  a 
machine  is  indicated,  iterations  proceed  from  the  feasible  revised  solu¬ 
tion  with  a  machine  added.  The  production  level  at  which  the  change 
occurs  becomes  the  upper  bound  of  the  production  range  for  the  base  solu¬ 
tion  and  the  lower  bound  for  the  feasible  range  of  the  revised  solution. 
Adding  any  number  of  machines  proceeds  in  the  same  fashion  since  adding 
8  hours  to  each  shift  has  the  same  effect  on  both  sides  of  equations 
(5.4)  and  (5.5),  regardless  of  the  current  level  of  machines. 
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All  aspects  of  the  decision  criteria  necessary  to  utilize  the 
dynamic  programming  formulation  have  now  been  discussed.  It  remains 
to  develop  a  computer  algorithm  which  incorporates  all  elements  of  the 
dynamic  and  linear  programming  solution  procedures.  This  development 
will  be  detailed  in  the  next  chapter. 


Chapter  VI 


DEVELOPMENT  OF  THE  COMPUTER  ALGORITHM  FOR  THE  MR  MODEL 

6.1  Introduction 

A  computer  algorithm  to  achieve  a  solution  to  the  dynamic  program¬ 
ming  formulation  of  the  machine  requirements  model  must  include  all  of 
the  decision  analysis  discussed  in  the  previous  chapter.  The  most  effec¬ 
tive  means  of  designing  such  an  algorithm  is  to  utilize  subroutines  to 
perform  those  aspects  of  the  procedure  which  recur  repetitively,  par¬ 
ticularly  the  linear  programming  iterations.  This  chapter  will  utilize 
a  discussion  and  information  flow  charts  to  explain  the  operation  of  the 
solution  procedure,  focusing  on  the  main  program  and  then  those  sub¬ 
routines  which  require  more  detailed  explanation.  Further  detail  is 
available  in  Appendix  B  where  the  complete  code  listing  is  included  along 
with  a  description  of  notation  used  in  the  FORTRAN  program. 

6.2  Main  Program 

The  main  program  of  the  computer  code  is  the  vehicle  which  coordin¬ 
ates  all  elements  of  the  solution  procedure.  Those  decision  points  which 
are  reached  within  a  subprogram  will  be  discussed  in  relation  to  that 
subroutine.  Figure  6.1  is  an  information  flow  chart  of  the  main  pro¬ 
gram.  Given  a  set  of  input  data  in  the  form  of  operating  rates,  per¬ 
cent  defectives,  production  costs  and  fixed  equipment  costs,  the  first 
step  in  the  solution  algorithm  Is  to  set  up  the  initial  simplex  tableau 
for  each  stage  of  the  production  process.  This  initial  tableau  will  be 
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Figure  6.1:  Information  flow  chart  of  main  program. 
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utilized  within  each  stage  as  it  is  reached  in  the  dynamic  programming 
formulation.  Following  this  initial  step,  the  most  cost  effective 
entering  variable  Is  determined  and  the  simplex  Iterations  within  each 
stage  begin.  Each  iteration  corresponds  to  a  distinct  production  range 
at  a  stage. 

Unless  there  is  a  range  split  due  to  a  change  in  cost  function  from 
the  previous  stage  or  a  machine  being  added,  each  change  In  production 
range  Is  marked  by  a  change  in  basis  requiring  the  formation  of  a  new 
tableau.  The  new  tableau  representing  revised  solution  information  is 
derived  in  the  algorithm  through  a  series  of  checks  and  matrix  manipula¬ 
tions.  Upon  determination  of  the  need  to  change  the  basis  to  attain  a 
higher  production  level,  the  PRELIM  subroutine  is  called.  This  routine 
determines  the  row  In  which  to  dual  simplex,  and  therefore  the  leaving 
variable,  and  finds  the  entering  variable  using  the  minimum  ratio  rule. 
Upon  establishing  the  new  basis,  the  INV  subroutine  is  called  to  invert 
the  basis  matrix  for  use  in  updating  the  tableau.  If  no  variable  is 
eligible  to  enter  the  basis,  indicating  all  of  the  feasible  rate  increases 
have  been  exhausted,  the  subroutine  to  add  a  machine,  MACHAD,  is  called 
and  the  resultant  tableau  is  then  utilized  for  further  Iterations. 

The  revised  basis  Inverse  can  now  be  used  in  conjunction  with  the 
initial  tableau  for  the  current  stage  to  obtain  the  updated  tableau. 

This  procedure  is  accomplished  in  the  UPDATE  subroutine  where  the  ini¬ 
tial  tableau  constraint  matrix  is  premultiplied  by  the  basis  inverse  to 
obtain  the  current  constraint  matrix  values.  The  right  hand  side,  ob¬ 
jective  function  and  Zj-Cj  values  are  also  updated  using  the  basis 
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Inverse  or  basic  variable  cost  values  where  appropriate.  The  updated 
tableau  now  contains  all  of  the  information  needed  to  determine  the 
production  range,  transition  function  values  and  applicable  cost  func¬ 
tion.  These  calculations  are  also  performed  within  the  UPDATE  routine. 

If  the  updating  of  a  tableau  is  necessary  within  a  dual  slmplexing  pro¬ 
cedure  designed  to  return  a  solution  to  feasibility,  a  RETURN  statement 
is  used  with  the  UPDATE  routine  to  bypass  the  calculations  not  associated 
directly  with  the  update  of  a  tableau.  Such  a  dual  simplex  procedure 
occurs  within  the  subroutines  which  check  for  adding  a  machine  (MACHCK) 
and  which  add  a  machine  when  rates  are  exhausted  (MACHAD). 

Each  time  a  new  solution  is  determined,  it  is  necessary  to  check  if 
adding  a  machine  may  be  economically  feasible.  If  this  check  is  positive, 
the  MACHCK  subroutine  is  called  to  make  this  analysis  as  explained  in 
Chapter  V.  Upon  this  determination,  the  necessary  updating  is  accom¬ 
plished  and  the  normal  procedure  continues  with  return  to  1  as  indicated 
in  the  flow  chart. 

When  the  upper  limit  of  production  is  reached  within  a  stage,  the 
solution  information  Is  compiled  and  printed  in  a  form  which  can  be  used 
to  develop  the  operational  data  which  optimizes  the  entire  production 
system.  Given  any  final  demand  value  for  production  within  the  upper 
limit  of  the  last  stage,  the  total  cost  of  production,  operating  rates 
at  each  stage,  number  of  machines  and  hours  of  operation  at  each  rate 
can  be  determined.  Two  of  the  subroutines  will  now  be  explained  In 
greater  detail. 
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6.3  Machine  Check  Subroutine 

The  flow  of  decisions  within  the  subroutine  MACHCK  for  determining 
the  economic  feasibility  of  adding  a  machine  is  depicted  in  Figure  6.2 
Since  the  procedure  of  the  analysis  Is  to  hypothetically  add  a  machine 
and  use  the  resultant  solution  for  further  analysis,  it  is  first  neces¬ 
sary  to  retain  all  pertinent  values  of  the  iteration  which  leads  to  the 
check.  This  Information  can  then  be  recovered  If  no  machine  is  added. 
When  a  machine  is  added  to  perform  the  analysis,  the  immediate  result 
is  that  the  solution  becomes  infeasible  for  the  range  comparison  needed. 
Therefore,  the  dual  simplex  procedure  must  be  used  to  drive  the  solution 
to  feasibility  at  which  time  the  ranges  will  be  compatible  for  analysis. 
This  dual  simplex  procedure  is  handled  by  the  RHSCK  subroutine  which  is 
used  solely  as  a  vehicle  to  determine  the  new  basis  while  iterating 
towards  feasibility.  Once  this  is  accomplished,  the  decision  criteria 
discussed  in  Section  5.4  is  utilized  to  decide  upon  the  addition  of  a 
machine  and  the  resultant  range  information.  Deciding  to  not  add  a 
machine  results  in  recovering  the  solution  tableau  which  initiated  the 
procedure  and  iterations  continue  from  there. 

6.4  Machine  Add  Subroutine 

The  MACHAD  subroutine  is  invoked  by  the  PRELIM  subroutine  when  it 
is  no  longer  feasible  to  bring  another  rate  into  the  basis.  Adding  a 
machine  makes  It  possible  to  achieve  higher  production  levels.  This  sub¬ 
routine  alters  the  right  hand  side  of  the  initial  tableau  for  the  current 
stage  to  reflect  the  Increase  in  production  hours  due  to  an  additional 
machine.  The  UPDATE,  RHSCK  and  INV  subroutines  are  then  used  to  drive 


Figure  6.2:  Information  flow  chart  of  the  MACHCK  subroutine 
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the  solution  to  feasibility  at  a  range  lower  bound  corresponding  to  the 
upper  bound  of  the  range  of  the  last  feasible  solution  without  the  addi 
tional  machine.  Continuity  of  production  ranges  is  maintained  and 
iterations  can  continue  from  this  point. 


Chapter  VII 


COMPUTATIONAL  RESULTS  AND  MODEL  VALIDATION 

7.1  Validation  of  the  Dynamic  Programming  Formulation 

In  order  to  be  assured  that  the  dynamic  programming  formula¬ 
tion  is  performing  all  of  the  functions  for  which  it  was  designed, 
it  is  necessary  to  test  the  model  against  a  valid  reference.  Although 
hand  calculations  can  be  used  to  verify  calculations  in  the  model, 
the  complexity  and  detail  of  the  solution  procedure  makes  such  a 
validation  prohibitive  for  all  but  the  most  simple  example.  Hand 
calculated  simplex  tableaus  were  used  to  test  some  elements  of  the 
model  but  the  most  meaningful  validation  results  from  using  a  commer¬ 
cial  integer  linear  programming  system  as  the  means  of  comparison  of 
results.  For  this  purpose,  the  Mathematical  Programming  System  (MPS) 
mixed  integer  linear  programming  package  (MISTIC)  was  utilized.  This 
system  is  available  through  the  Virginia  Tech  Computing  Center. 

A  four  stage  manufacturing  process  with  input  data  as  shown  in 
Table  7.1  was  used  as  the  validation  example.  The  procedure  of  using 
the  results  of  the  dynamic  programming  solution  algorithm  is  shown  in 
Table  7.2.  This  table  depicts  the  production  range  of  each  stage  and 
the  resulting  operational  information  for  optimization  of  producing 
to  a  final  demand  of  500  units.  The  procedure  starts  with  the  fourth 
stage  since  the  final  production  Is  known  and  the  decisions  can  be 
determined  through  recursive  movement  through  the  stages.  The  complete 
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Table  7.1  Input  Data  for  4  Stage  Production  Process 


Stage  Production  Data 


1 

Rates  (Units/hour) 

25 

30 

35 

40 

45 

50 

%  Defectives 

.041 

.042 

.045 

.050 

.051 

.055 

1st  Shift  Cost 
($/hour) 

10.0 

12.0 

14.0 

17.2 

19.8 

25.0 

2nd  Shift  Cost 
($/hour) 

13.0 

13.8 

16.1 

19.2 

22.9 

30.0 

2 

Rates 

10 

15 

20 

25 

30 

%  Defectives 

.094 

.095 

.101 

.105 

.110 

1st  Shift  Cost 
($/hour) 

7.0 

7.95 

9.00 

10.75 

13.2 

2nd  Shift  Cost 
($/hour) 

10.00 

10.95 

13.00 

14.00 

17.40 

3 

Rates 

40 

45 

50 

55 

60 

%  Defectives 

.045 

.047 

.050 

.051 

.053 

1st  Shift  Cost 
($/hour) 

5.2 

5.85 

8.0 

9.35 

10.8 

2nd  Shift  Cost 
($/hour) 

7.2 

7.65 

10.0 

11.0 

13.8 

4 

Rates 

5 

10 

15 

20 

25 

1  Defectives 

.057 

.060 

.063 

.067 

.072 

1st  Shift  Cost 
($/hour) 

10.0 

15.0 

19.95 

30.0 

40.0 

2nd  Shift  Cost 
($/hour) 

15.0 

20.0 

30.0 

35.0 

50.0 

Fixed  Cost  of  Equipment  Per  Day: 

Stage  Cost  (S/day) 

1  80.0 

2  60.0 

3  40.0 

4  70.0 


59 


solution  results  for  the  four  stage  example  are  shown  in  Appendix  C 
and  these  can  be  used  to  find  the  optimal  operational  data  for  any 
final  demand  up  to  753  units. 

To  test  the  results  obtained  from  the  dynamic  programing  model 
and  to  assure  that  all  routines  of  the  model  are  exercised,  two  valida¬ 
tion  computer  runs  were  made  using  the  MISTIC  procedure.  One  test  was 
for  a  final  demand  of  500  units  and  the  other  for  a  final  demand  of  650 
units.  Table  7.3  depicts  the  results  of  the  comparison  of  the  two 
models  for  each  of  the  demand  values.  These  results  verify  that  both 
procedures  reach  virtually  identical  results  for  a  specific  final 
demand  when  maintaining  an  integer  number  of  machines.  Of  particular 
note  is  the  difference  in  execution  time  to  attain  the  same  solution 
using  the  different  procedures.  Not  only  is  the  dynamic  programming 
formulation  much  faster,  using  an  IBM  370/158  computer,  but  it  only  has 
to  be  executed  once  to  obtain  solution  information  for  any  projected 
final  demand  up  to  the  upper  limit  planned  for.  The  integer  linear  pro¬ 
gramming  procedure  must  be  executed  for  each  distinct  final  demand  value. 

7.2  Computational  Efficiency  of  the  Model 

Several  tests  of  the  dynamic  programming  model  formulation  were 
conducted  to  determine  the  computational  effect  as  related  to  the  MISTIC 
procedure.  In  addition  to  the  four  stage  process  already  examined,  the 
number  of  operating  rates  for  the  four  stage  example  were  doubled  to 
test  the  effect  of  increased  number  of  variables  on  the  decision  pro¬ 
cess.  To  analyze  the  effect  of  an  increased  number  of  stages  in  the 
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manufacturing  process,  a  six  stage  and  a  five  stage  example  were  for¬ 
mulated  and  tested.  All  of  the  examples  were  solved  using  both  the  DP 
procedure  and  MISTIC.  The  results  of  these  test  problems  are  shown  in 
Table  7.4  and  the  production  costs  are  identified  for  each  method  to 
further  support  the  credibility  of  the  DP  model. 

The  most  significant  aspect  of  the  comparative  results  is  the 
major  decrease  in  computation  time  for  the  DP  algorithm  as  compared 
to  the  integer  programming  method.  As  the  number  of  stages  increases, 
the  time  to  attain  a  solution  with  the  MISTIC  procedure  rapidly  in¬ 
creases  while  the  DP  model  leads  to  a  more  linear  increase.  The  model¬ 
ling  of  even  higher  numbers  of  stages  should  realize  an  even  greater 
discrepancy  in  execution  times.  Such  a  result  indicates  that  the 
MISTIC  code  or,  in  fact,  any  branch  and  bound  mixed  integer  procedure, 
would  be  prohibitive  for  a  complex  manufacturing  process. 

Two  other  areas  for  consideration  in  comparing  the  two  solution 
procedures  are  input  data  requirements  and  the  computer  core  storage 
needed.  The  data  requirements  for  a  programming  package  such  as  MPS 
are  considerably  more  complex  than  that  of  the  DP  model  as  demonstrated 
In  Table  7.1.  MPS  is  limited  to  two  data  iter,<  per  card  and  the  task 
of  preparing  the  data  cards  for  a  large  linear  programming  constraint 
matrix  is  monumental.  The  DP  formulation  contains  routines  to  form  the 
linear  programming  tableaus  from  simple  production  data  and  much  prepa¬ 
ration  time  and  sources  for  error  are  eliminated.  Table  7.4  also  Indi¬ 
cates  the  core  requirements  for  the  computer  runs  of  the  example  pro¬ 
blems.  Much  less  core  is  required  in  the  DP  algorithm  resulting  in 
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application  to  computer  systems  with  smaller  memory  capacity  and  less 
cost  of  computer  operation. 

The  significant  results  obtained  through  solving  the  machine  re¬ 
quirements  problem  with  dynamic  programming  offers  a  base  upon  which 
to  build  further  research  In  manufacturing  systems  optimization.  Ex¬ 
tension  of  the  model  Into  additional  areas  will  be  discussed  in  the  next 
chapter. 
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Chapter  VIII 


EXTENSIONS 

8,1  Extensions  Related  to  Dynamic  Programming  Formulation 

The  dynamic  prograitmi ng  formulation  of  the  machine  requirements 
model  lends  itself  to  several  variations  on  the  model  depicted  in 
Chapter  IV.  Any  manufacturing  process  which  can  be  handled  in  distinct 
stages  and  which  has  similar  production  and  cost  data  can  be  modelled. 
This  feature  allows  extension  of  the  model  to  non-serial  systems  such 
as  converging  and  diverging  branch  and  feedforward  and  feedback  systems. 
In  essence,  non-serial  systems  can  be  decomposed  into  a  series  of  serial 
systems  for  solution  purposes.  The  extension  of  the  model  to  these  types 
of  systems  would  add  complexity  to  the  absorption  of  previous  stage  pro¬ 
duction  cost  as  well  as  to  the  transition  function.  The  basic  solution 
procedure  is  capable  of  handling  such  complexities  within  the  linear 
programming  structure  at  each  stage. 

The  feasibility  of  extending  the  DP  formulation  to  non-serial  sys¬ 
tems  allows  the  procedure  to  be  applied  to  complex  manufacturing  systems 
having  a  variety  of  operating  modes.  Such  production  situations  as  manu¬ 
facturing  twoormore  products  simultaneously  for  all  or  a  portion  of  a 
production  line,  rework  of  defective  Items,  and  the  converging  or  diverg¬ 
ing  of  machining  processes  in  a  system  can  all  be  treated  within  the 
framework  of  non-serial  systems.  More  detail  on  the  structure  and  solu¬ 
tion  of  such  systems  is  available  in  Nemhauser  [24]. 
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An  example  of  a  non-serial  system  which  is  particularly  amenable 
to  a  dynamic  programming  solution  procedure  is  the  converging  branch  sys¬ 
tem  as  diagrammed  in  Figure  8.1.  This  system  consists  of  three  serial 
flow  machining  systems  which  converge  into  a  final  assembly  point.  Since 
quantity  input  to  the  last  stage  must  be  the  output  from  the  previous 
stage,  or  stages,  in  the  assembly  situation  as  shown,  the  output  from 
the  three  converging  branches  must  be  equal.  This  simplifies  the  pro¬ 
blem  considerably  as  each  converging  branch  can  be  solved  individually 
with  the  output  of  the  branch  required  to  be  the  input  needed  by  the 
last  stage  for  a  specified  final  demand  value.  The  total  cost  of  the 
system  through  the  assembly  process  would  be  composed  of  the  sum  of  all 
of  branches  and  the  cost  of  final  assembly. 

The  development  of  the  dynamic  programming  model  has  focused  on  a 
manufacturing  system  with  the  stages  limited  to  machining  centers.  Such 
a  concept  proceeds  directly  from  the  historical  perspective  of  the  pro¬ 
blem  and  the  practical  consideration  of  limiting  the  scope  of  the  problem 
to  the  most  relevant  aspects.  This  approach  is  also  consistent  with  the 
incorporation  of  the  machine  parameters  problem  into  the  machine  require¬ 
ments  problem  structure  where  all  machining  processes  are  modelled  as 
having  discrete  operating  rates  of  operation  and  associated  costs  can 
be  built  Into  the  optimization  procedure. 

A  particular  example  of  modelling  a  non-machining  process  would  be 
material  handling  between  machining  stages.  In  the  design  of  a  manu¬ 
facturing  process,  the  type  and  number  of  material  handling  equipment 
may  be  critical.  Given  a  set  of  alternative  equipment  with  discrete 


Figure  8.1:  Converging  branch  manufacturing  system. 
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operating  rates,  these  rates  could  be  used  to  develop  the  decision  vari¬ 
ables  for  a  material  handling  stage  between  the  machining  stages.  The 
operating  cost  per  unit  time  must  be  known  as  well  as  the  fixed  equipment 
cost.  The  development  of  such  an  integrated  concept  would  take  the  op¬ 
timization  of  the  manufacturing  process  a  step  further  and  limit  the 
variability  associated  with  actual  implementation  of  the  model  due  to 
restrictions  imposed  by  materials  handling  equipment  selection. 

8.2  Extensions  Related  to  Model  Structure 

In  addition  to  extensions  related  to  the  DP  solution  procedure,  the 
structure  of  the  machine  requirements  model  is  such  that  a  variety  of 
situations  can  be  incorporated.  The  model  developed  in  Chapter  VI 
utilizes  a  two  shift  production  system  but  the  time  system  employed  in 
the  model  is  highly  flexible.  The  use  of  such  modes  as  3  shifts,  2 
shifts  with  overtime  and  limitations  to  overtime  can  be  readily  achieved 
through  variations  in  the  constraints  or  the  addition  of  new  constraints. 
These  constraint  changes  then  lead  to  revising  the  set  of  decision  vari¬ 
ables  In  much  the  same  way  as  increasing  the  number  of  rates  or  stages 
affected  the  test  problem  discussed  in  Chapter  VII. 

Two  elements  of  the  machine  requirements  model  developed  by  Davis 
and  Miller  [7]  can  be  included  in  the  dynamic  programming  formulation. 

The  Davis  and  Miller  model  is  Identical  to  the  model  developed  in  this 
research  except  for  their  Inclusion  of  in-process  inventories.  Also, 
their  model  uses  units  processed  as  the  decision  variables  in  each  stage 
rather  than  the  hours  used  for  production,  but  there  is  a  direct 
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relationship  between  these  variables.  Through  analysis  of  the  Davis  and 
Miller  model  it  can  be  seen  that  the  effect  of  in-process  inventories 
ultimately  appears  only  in  the  objective  function  where  storage  cost  can 
directly  affect  the  cost  Incurred  from  a  certain  operating  rate.  There¬ 
fore,  the  DP  formulation  can  be  altered  to  include  in-process  inventories 
through  modification  of  the  cost  data.  In  this  way,  the  model  can  be 
designed  to  accommodate  the  cost  of  storage  between  stages  as  it  affects 
the  cost  of  the  entire  manufacturing  system. 

The  Davis  and  Miller  model  is  also  designed  to  handle  discretely 
distributed  production  demand  through  the  use  of  probability  data  for 
distinct  demand  values  and  cost  data  for  shortage  and  overage  produc¬ 
tion.  Such  a  feature  can  be  readily  added  to  the  DP  model  if  the  proba¬ 
bility  information  is  available,  since  the  calculations  occur  outside  of 
the  main  optimization  process.  The  DP  formulation  actually  makes  this 
feature  more  readily  applicable  since  all  of  the  optimization  results 
are  available  for  an  entire  range  of  production.  The  optimum  cost  for 
producing  to  each  specific  demand  value  can  be  easily  calculated  as 
shown  in  Table  7.2  and  then  the  probability  and  cost  data  can  be  applied 
to  determine  the  actual  expected  cost.  Recourse  to  an  integer  formula¬ 
tion  such  as  MISTIC  requires  a  complete  solution  procedure  for  each  dis¬ 
tinct  demand  value  for  which  probability  information  is  available. 

Although  the  main  approach  in  this  research  has  been  to  consider 
the  optimization  of  a  machining  system  from  the  design  standpoint,  the 
availability  of  solution  data  for  an  entire  range  of  production  readily 
allows  extension  to  a  production  scheduling  procedure.  The  fixed  cost 


70 


of  a  machine  in  the  model  formulation  includes  Investment  cost  since  the 
number  of  machines  is  being  planned.  For  an  optimization  procedure  of 
a  production  line  already  installed,  the  investment  portion  of  the 
fixed  daily  cost  is  removed.  An  upper  bound  must  also  be  placed  on  the 
number  of  machines  allowed  to  be  allocated  to  a  machining  center.  The 
production  time  decision  variables  can  readily  be  transformed  to  produc¬ 
tion  quantity  if  desired  by  production  personnel. 


Chapter  IX 


SUMMARY  AND  CONCLUSIONS 

In  summary,  it  has  been  shown  that  it  is  possible  to  discretize 
the  machining  parameter  problem  to  enable  any  type  of  machining 
equipment  to  be  incorporated  into  the  machine  requirements  planning 
problem.  The  machining  parameter  algorithm  is  more  efficient  than 
any  other  developed  in  the  literature  and  is  the  only  one  which  results 
in  discrete  operational  data  as  shown  in  Table  3.2.  The  dynamic 
programming  formulation  of  the  machine  requirements  model  has  been 
shown  to  lead  to  a  highly  efficient  and  flexible  solution  procedure. 

The  computational  efficiency  as  depicted  in  Table  7.4,  results  in  the 
ability  to  optimize  complex  multistage  systems  without  incurring 
prohibitive  solution  costs  and  core  storage  requirements. 

In  many  complex  systems,  branch  and  bound  integer  programming 
algorithms  may  prove  to  be  incapable  of  reaching  a  solution  in  a 
feasible  amount  of  time.  Thus  a  DP  approach  is  the  only  recourse 
if  an  optimal  solution  is  to  be  obtained.  The  fact  that  solution 
information  for  an  entire  range  of  production  demand  is  available 
also  contributes  greatly  to  reducing  computational  time  when  optimizing 
for  a  set  of  final  demand  values.  The  feasibility  of  the  dynamic 
programming  approach  as  demonstrated  by  the  research  can  lead  to  a 
much  greater  use  of  optimization  techniques  in  manufacturing  processes. 

The  development  of  the  discrete  approach  to  the  machining 
parameter  problem  is  the  key  aspect  of  the  research  which  allows  the 
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DP  model  to  be  applied  to  either  AC  or  DC  machining  systems.  In 
addition,  machining  systems  with  any  conti nous  power  source,  such  as 
hydraulic  or  AC  systems  with  a  continuously  variable  drive,  can  be 
discretized  with  the  MP  algorithm.  Through  the  use  of  a  basic, 
practical  assumption,  discrete  operating  rate  and  cost  information 
is  attainable  for  any  machine  center  with  continuous  drive  equipment. 
This  data  can  then  be  incorporated  into  the  stages  of  the  MR  model 
resulting  in  a  direct  interface  between  the  two  optimization  procedures. 
The  importance  of  the  discrete  nature  of  the  MP  formulation  becomes 
clear  with  the  realization  that  a  high  percentage  of  machinery  is  of 
a  continuous  nature.  Although  the  machining  parameters  model  can 
optimize  the  parameters  for  a  particular  process,  this  aspect  of  the 
procedure  is  unnecessary  for  the  interface  with  the  MR  model.  All 
feasible  operating  rate  data  is  needed  to  optimize  the  total  system. 
Attempting  to  limit  the  data  to  only  rates  which  can  be  expected  to 
be  utilized  would  not  significantly  reduce  the  complexity  of  the 
solution  procedure.  This  is  supported  by  the  results  of  the  test 
to  analyze  the  effect  of  increasing  the  number  of  rates.  Only  a 
slight  increase  in  solution  time  was  exhibited  for  the  DP  solution. 

A  distinct  advantage  of  the  dynamic  programming  approach  to  the 
MR  problem  is  the  variety  of  machining  systems  which  can  be  effect¬ 
ively  optimized.  The  extension  of  the  model  to  non-serial  systems 
without  greatly  increasing  the  complexity  of  the  solution  procedure 
allows  optimization  techniques  to  be  applied  to  virtually  all  types  of 
machining  systems.  The  flexibility  of  the  model  is  also  demonstrated 
by  the  ease  with  which  different  aspects  such  as  changes  in  time 
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periods,  demand  probability  data,  and  cost  of  storage  between  stages 
can  be  incorporated. 

An  integrated  manufacturing  system  optimization  model  as  developed 
in  this  research  contributes  significantly  to  machining  technology 
state  of  the  art.  The  systems  approach  using  dynamic  programming 
should  lead  to  modeling  additional  aspects  of  the  manufacturing 
process.  With  the  development  of  efficient  models  which  expand  into 
all  realms  of  manufacturing  industry,  the  use  of  computers  and  auto¬ 
mation  can  reach  its  full  potential  for  improving  industrial  product¬ 
ivity. 
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APPENDIX  C 

OPERATIONAL  DATA  RESULTS  FOR  MACHINE  REQUIREMENTS 
PLANNING  MODEL  USING  4  STAGE  EXAMPLE 
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NOTATION  FOR  THE  INTERPRETATION  OF  HEADINGS  FOR  RESULTS 


-  Number  of  the  machining  center  of  the  system 

-  Lower  Bound  of  the  production  range  for  the  current  solution 

-  Upper  bound  of  the  production  range 

-  Number  of  machines  at  the  machining  center 

-  Fixed  value  of  the  transition  function  which  specifies  the 
input  to  the  current  stage  and  the  output  from  the  previous 
stage 

-  Variable  value  of  the  transition  function  which  is  multiplied 
by  the  production  output  required 

-  Fixed  portion  of  the  cost  function  for  producing  a  specified 
output  through  the  current  stage 

-  Variable  portion  of  the  cost  function  which  must  be  multi¬ 
plied  by  the  production  output  required 

-  Index  number  of  the  basic  variable  which  is  analogous  to  the 
rate  of  operation.  The  variables  are  indexed  through  twice 
the  number  of  operating  rates  to  account  for  both  first  and 
second  shifts.  Second  shift  rates  are  Indicated  by  index 
numbers  which  begin  one  greater  than  the  number  of  operating 
rates  for  the  machine. 

-  The  operating  rate  of  the  machine 

-  Fixed  portion  of  the  amount  of  time  in  hours  which  the 
machine  should  be  operated  at  the  specified  rate 

-  Variable  amount  of  time  for  which  the  machine  Is  to  be 
operated.  This  value  must  be  multiplied  by  the  output 
needed  from  the  stage  and  added  to  FIXTIME  to  obtain  the 
total  time  of  operation.  Total  units  processed  at  the  stage 
can  also  be  obtained  by  multiplying  the  time  by  the  operat¬ 
ing  rate  and  summing  for  each  stage. 


0000*0  0**0tt  t*0O*l  00*0  l  00*0001  00*000 
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A  GENERALIZED  MACHINE  REQUIREMENTS  PLANNING  ALGORITHM 
FOR  SERIAL  FLOW  MACHINING  SYSTEMS 


by 

Glenn  M.  Hayes,  Jr. 

(ABSTRACT) 

The  machine  requirements  planning  problem  is  one  which  is  appli¬ 
cable  throughout  the  manufacturing  industry.  The  use  of  automation 
and  computer  technology  makes  the  use  of  machine  requirements  modelling 
particularly  attractive,  and  the  vast  capital  investment  in  machining 
equipment  offers  significant  opportunity  for  savings  through  optimi¬ 
zation.  The  machine  requirements  problem  must  be  modelled  as  a  mixed 
integer  linear  program.  However,  a  dynamic  programming  solution  pro¬ 
cedure  with  linear  programming  postoptimality  techniques  at  each  stage 
of  the  machining  process  is  utilized  to  provide  an  efficient,  flexible 
algorithm.  The  machine  requirements  model  is  capable  of  optimizing 
any  type  of  machining  system,  whether  of  a  discrete  or  continuous 
operating  mode.  Through  discretizing  the  machine  parameters  problem, 
by  limiting  tool  changes  to  between  passes  of  a  tool,  any  continuous 
mode  system  Is  approximated  as  discrete  for  optimization  purposes  and 
a  direct  Interface  with  the  machine  requirements  planning  model  is  pro¬ 
vided. 

The  dynamic  programming  solution  procedure  is  compared  with  a  mixed 
Integer  procedure.  The  DP  formulation  is  not  only  more  efficient  in 


l 
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both  time  and  core  but  provides  sensitivity  Information  and  offers 
broad  spectrum  of  further  application  Into  more  complex  aspects  of 
manufacturing  systems. 


